Sensitivity to neutrino mixing parameters with atmospheric neutrinos 



Abhijit Samanta * 

Harish- Chandra Research Institute, Chhatnag Road, Jhusi, Allahabad 211 019, India 

(Dated: November 19, 2009) 

We have analyzed the atmospheric neutrino data to study the octant of #23 and the precision of the 
oscillation parameters for a large Iron CALorimeter (ICAL) detector. The ICAL being a tracking 
detector has the ability to measure the energy and the direction of the muon with high resolution. 
From bending of the track in magnetic field it can also distinguish its charge. We have generated 
events by Nuance and then considered only the muons (directly measurable quantities) produced in 
charge current interactions in our analysis. This encounters the main problem of wide resolutions 
of energy and baseline. The energy-angle correlated two dimensional resolution functions are used 
to migrate the energy and the zenith angle of the neutrino to those of the muon. A new type of 
binning has been introduced to get better reflection of the oscillation pattern in chi-square analysis. 
Then the marginalization of the x 2 over all parameters has been carried out for neutrinos and anti- 
neutrinos separately. We find that the measurement of #13 is possible at a significant precision with 
atmospheric neutrinos. The precisions of Am^ and sin 2 #23 are found ~ 8% and 38%, respectively, 
at 90% CL. The discrimination of the octant as well as the deviation from maximal mixing of 
atmospheric neutrinos are also possible for some combinations of (023, #13 )■ We also discuss the 
impact of the events at near horizon on the precision studies. 

PACS numbers: 14.60.Pq 



I. INTRODUCTION 



Recent discovery of the neutrino mass has opened up 
a new window into physics beyond the standard model. 
Aside from this fact, the two surprising sets of results [lj], 
i) the extremely small masses of neutrinos (very different 
from quark sector) and ii) a dramatically different mix- 
ing pattern from quarks, indicate a new direction of this 
field. The first one may be the hint of a new symme- 
try such as B — L at high scales so that one can use a 
mechanism like seesaw to resolve the puzzle of the small- 
ness of the masses. On the other hand, the second one 
poses a much more challenging problem. One can ex- 
pect a new symmetry for leptons as well as for quarks to 
solve this problem. Currently, there are many theoretical 
ideas. For example, the \i — t symmetry [2| is invoked 
to explain the maximal mixing. However, if this \i — t 
symmetry exists, it leads to 8c p = and 6*13 = 0. It 
should be broken if there appears a nonzero #13 and CP 
violation. In that case the octant (the sign of 6*23 — 45°) 
and the nonzero value of #13 emerges other new possibil- 
ities. It is also expected that neutrino theories may have 
implications on the very fascinating fields like observed 
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matter-antimatter asymmetry of the universe, grand uni- 
fication, supersymmetry, extra dimensions, etc[3|. 

Active endeavors are under way to launch the era of 
precision experiments with a thrust to uncover the un- 
derlying principle that gives neutrino masses and their 
mixing. This is one of the most promising ways to ex- 
plore physics beyond the standard model. In the stan- 
dard oscillation picture there are six parameters. The 
present la, 2a and 3a confidence level (CL) ranges from 
global 3v oscillation analysis (2008) [l[ are very exciting 
(see table Recently, new bounds, #13 = —0.07^°;^ 
and the asymmetry 6» 23 - tt/4 = 0.03±g;?§ at 90% CL 
have been shown in 0, Q from an analysis considering 
all present neutrino data. The CP-violating phase 8c p 
is still unconstrained. 

Despite of these spectacular achievements, a lot of 
things are still missing. Tremendous efforts are underway 
to determine the mass ordering (sign of Am| 2 ), the values 
of 0i3 and Sep, and to discriminate the octant degeneracy 
of #23 in future experiments. We define Am| 2 = m|— mi- 
There are many ongoing and planned experiments: UNO 
0, T2K 0, NOvA [jj, Hyper-Kamiokande 0, INO [l0| 
and many others. The main characteristic feature of a 
magnetized Iron CALorimeter (ICAL) detector proposed 
at India-based Neutrino Observatory (INO) is that it has 
the capability to detect and separately, which mea- 
sures directly the matter effect. 



2 



Parameter Am^/ICT 5 eV 2 sin 2 12 sin 2 6 13 sin 2 6 23 | Am§i |/10~ 3 eV 2 

Best fit 7^67 0.312 0.016 0.466 2^39 

lcr range 7.48 - 7.83 0.294 - 0.331 0.006 - 0.026 0.408 - 0.539 2.31 - 2.50 

2a range 7.31 - 8.01 0.278 - 0.352 < 0.036 0.366 - 0.602 2.19 - 2.66 

3<t range 7.14 - 8.19 0.263 - 0.375 < 0.046 0.331 - 0.644 2.06 - 2.81 



TABLE I: Global 3u oscillation analysis (2008) 



Unlike a fixed baseline neutrino beam experiment, the 
atmospheric neutrino flux covers a wide range of baseline 
(a few km - 12900 km) and energy (sub GeV - a few 
hundred GeV). On the other hand, it is not known well 
and there are huge uncertainties in its estimation. It is 
also a very rapidly falling function of energy. So, the 
extraction of the results from the experimental data is 
very complicated. 

The deviation from maximal mixing and the discrimi- 
nation of octant degeneracy of #23 have been studied in 
[ill [l2| with atmospheric neutrinos for a large magne- 
tized ICAL detector. However, the results have been ob- 
tained without marginalization and assuming the Gaus- 
sian resolution functions of fixed widths for whole range 
of energy and zenith angle. The energy range for the 
atmospheric neutrinos is very wide. The resolutions are 
changed significantly over its range and are very different 
for neutrinos and anti-neutrinos. Moreover, the energy 
resolutions appear to be non-Gaussian due to some un- 
measurable product particles like neutral hadrons in neu- 
trino interaction even if one considers all visible hadrons. 

For a given neutrino energy and direction, there is a 
distribution in the reconstructed energy and direction. 
Again, a particular reconstructed energy and direction 
can come from a wide range of true neutrino energy and 
direction. So, it is not possible to convert a distribu- 
tion in reconstructed energy and direction obtained from 
an experiment to a distribution in actual neutrino en- 
ergy and direction. This restricts the binning of the data 
for chi-square analysis only in experimentally measured 
energy and direction. On the other hand, the actual reso- 
lution functions have no regular pattern and significantly 
deviate from the Gaussian nature even if we consider the 
visible hadrons. Again, the width changes with neutrino 
energy. For a simplistic analysis, if one considers a Gaus- 
sian resolution with a width that gives equal space under 
the surface of resolution function, the correct theoreti- 
cal data smearing this approximated Gaussian resolution 
function can not be obtained for chi-square analysis. As 
a consequence, the best-fits and the contours of oscilla- 



tion parameters will differ largely from the true values. 
In literature, there are many analyses where both the 
theoretical as well as the experimental data are obtained 
by smearing the Gaussian resolution functions. For an 
example, see ref. 12j. However, the result changes very 
rapidly with change of the width of the resolution. So, 
realistic estimation of the capability of an experiment can 
be done only by an analysis with experimentally measur- 
able quantities and exact resolution functions. 

Till now, the precision studies with atmospheric neu- 
trinos have mainly carried out for water Cherenkov de- 
tector, a non-magnetized detector. It is very important 
to see the capability of a large magnetized detector. We 
have studied the neutrino oscillation considering neutri- 
nos and anti-neutrinos separately in the chi-square anal- 
ysis. Here, we consider the muons (directly measurable 
quantities at ICAL) produced by the charge current in- 
teractions. We generate events by Nuance- v3 13|. The 
two dimensional energy-angle correlated resolution func- 
tions are used to migrate the energy and the zenith angle 
of the neutrino to the energy and the zenith angle of the 
muon. 

The above method has been introduced in 1141 and 
later used in 15 1 . The goal of the previous work [15j was 



only to compare the allowed parameter space of oscilla- 
tion parameters obtained from different types of binning. 
The considered systematic uncertainties were very much 
different from the present systematic uncertainties. The 
purpose of this work is to study the following. 

We consider whole data set in previous studies. But in 
reality, the horizontal events cannot be detected when the 
iron slabs are stacked horizontally. In this paper, we have 
studied the impact of these events in determining the 
precision of the parameters with and without considering 
a rejection criteria for the horizontal events. This is very 
crucial to determine whether horizontal stacking of iron 
plates is better than the vertical stacking or not. 

As discussed in 15[ , the binning of the data neither in 



log E — cos ^zenith nor in log E — log L is the optimum. In 
this paper, we have optimized the binning in L. These 
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are equal binned grids in log E — L 0A plane, which can 
capture the oscillation behavior for all L and E in a better 
way in the chi-square analysis. Again, the number of 
bins in both axes need optimization between resolutions 
and statistics. However, it should be noted here that if 
the statistics is huge for whole range of E and L, one 
can solve this problem by making the bin size very small 
and then the type of binning will not play any crucial 
role. However, the type of binning is very crucial when 
the analyses is in experimentally measurable energy and 
directions. Here the statistics over measured energy and 
direction is redistributed notably from the true neutrino 
energy and direction. 

Finally, we have made a detailed study on the sensi- 
tivity of a magnetized ICAL detector in determining the 
precision of Am^ and #13 as well as in discriminating 
the octant ambiguity of 6*23. We find the sensitivities of 
the parameters in two dimensional parameter space after 
marginalization over whole allowed ranges of the param- 
eters. The absolute bounds of each parameter are also 
studied. 



II. ATMOSPHERIC NEUTRINO FLUX AND 
EVENTS 

The atmospheric neutrinos are produced by the in- 
teractions of the cosmic rays mainly with nucleuses of 
molecules in the earth's atmosphere. The knowledge of 
primary spectrum of the cosmic rays has been improved 
from the observations by BESs[li| and AMS[17]. How- 
ever, large regions of parameter space have not been ex- 
plored and they are interpolated or extrapolated from 
the measured flux. The difficulties and the uncertain- 
ties in the calculation of the neutrino flux depend on 
the neutrino energy. The low energy fluxes have been 
known quite well. The cosmic ray fluxes (< 10 GeV) 
are modulated by the solar activity and the geomagnetic 
field through a rigidity (momentum/charge) cutoff. At 
the higher neutrino energy (> 100 GeV), the solar ac- 
tivity and the rigidity cutoff are irrelevant [3| . There is 
10% agreement among the calculations for neutrino en- 
ergy below 10 GeV because different hadronic interaction 
models are used in the calculations and because the un- 
certainty in the cosmic ray flux measurement is 5% for 
the cosmic ray energy below 100 GeV [ijj]. In our sim- 
ulation, we have used a typical Honda flux calculated in 
3-dimensional scheme 



rial are simulated using the Monte Carlo model Nuance 
(version- 3 ) [3| . Here, the charged current (CC) and neu- 
tral current (NC) interactions are considered for (quasi- 
elastic, resonance, coherent, diffractive, and deep inelas- 
tic scattering processes. 



III. OSCILLATION OF ATMOSPHERIC 
NEUTRINOS 

The present atmospheric neutrino data are well ex- 
plained by two flavor oscillation 19|, l20|. However, one 
expects a considerable — > v e oscillation of atmospheric 
neutrinos in 3-flavor framework if $13 is nonzero. To un- 
derstand the analytical solution one may adopt the so 
called "one mass scale dominance" (OMSD) frame work: 
lAm^l << |mj — m 2 2 |- Then the oscillation probabili- 
ties can be expressed as: 
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These oscillation probabilities are derived for vacuum. 
Since the oscillation involves electron neutrino, it will be 



modulated by the matter effect 2lj, l22| . Then 
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Here, E, L and Amf^ are in GeV, km and eV 2 , respec- 
tively. 
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The interactions of neutrinos with the detector mate- 



The mass squared difference (Amf^)" 1 and mixing an- 
gle sin 2 26*^3 in matter are related to their vacuum values 
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FIG. 1: The oscillogram of — » oscillation probability in E — cos ZO nith plane for 823 = 40° (left column) and 50° (right column) 
with 81s — 5° (lower row) and 7.5° (upper row). We choose Am^ = —2.5 x 10 -3 eV 2 and 8cp = 0. 



by 



(Am^)™ = ^{(Am 2 31 ) cos26>i 3 - A) 2 + {{Am 2 31 ) sm29 13 ) 

(Am^)sin 29 13 



[3] and [4] that a resonance in above oscillation probabil- 
ities will occur for neutrinos (anti-neutrinos) with NH 
(IH) when 



sin29 



13 
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=(4) 



sin 2 26™3 -> 1 or, A = &m\ x cos 26»i 3 . 



where, A = 2^/2GfN c E, Gf is the Fermi constant, 
N e is the electron density of the medium and E is neu- 
trino energy 23[. The matter potential term A has 
the same absolute value, but opposite sign for neutrino 
and anti-neutrino. The superscript 'm' denotes effective 
parameters in matter. Due to this matter effect, the 
Mikheyev-Smirnov- Wolfenstein (MSW) resonance occurs 
in P(Ufj, — > v e ) or ¥(v e — * v^). It happens for Normal Hi- 
erarchy (NH) with neutrinos and for Inverted Hierarchy 
(IH) with anti-neutrinos. It can be understood from Eq. 



Then the resonance energy can be expressed as 



E = 



1 



2 x 0.76 x 10- 4 K, 



eV 2 



13 



gm/cc 



(5) 



.(6) 



The resonance energy corresponding to a baseline can be 
seen in 



14] 



The oscillogram of muon survival probability is demon- 
strated in Fig. CDfor 9 13 = 5° and 7.5° with 6> 23 = 40° and 
50°, respectively. Here, we show the resonance ranges 
for the neutrinos passing through the core of the earth 
(with Ek3-6 GeV) and the mantle of the earth (with 
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E « 5 - 10 GeV). We also see a difference for 6» 23 = 40° 
and 50° due to the sin #23 term (Eq. [3]) , which domi- 
nates over the other terms due to the matter effect. 

IV. BINNING OF THE EVENTS 

For binning of the data in E, we need to consider the 
following facts. I) The atmospheric neutrino flux falls 
very rapidly with increase in energy. II) Again, the wide 
resolutions of E and L between true neutrinos and recon- 
structed neutrinos smear the oscillation effect to a signifi- 
cant extent. The wide resolutions arise mainly due to the 
interaction kinematics. This huge uncertainty in recon- 
structed neutrino momentum is due to the un-observable 
product particles and slightly due to the un-measurable 
momentum of recoiled nucleus when E ^ 1 GeV. For this 
reason, the energy resolutions deviate largely from the 
Gaussian nature. These are strongly neutrino energy de- 
pendent. At low energy (E ^ 1.5 GeV) the quasi-elastic 
process dominates and the muon carries almost whole en- 
ergy of the neutrino. The energy resolution is very good 
here. With increase in energy, the width of the resolution 
increases significantly as the deep inelastic event domi- 
nates as well as the flux also falls very rapidly. This is 
one of the main problems in the atmospheric neutrino 
experiments. Ill) There is also an important characteris- 
tics of the oscillation probability when both L and E are 
varied simultaneously. We explain it here for «-> 
oscillation in vacuum, which is a sinusoidal function of 
L/E. If we plot it in L — E plane (see Fig. |2J), it is 
seen that the distance between two consecutive peaks of 
oscillation in E for a fixed L increases very rapidly with 
E. These three points suggest increase in bin size with 
increase in E. We choose equal bin size in log-E 1 . 

Again, the distance between two consecutive peaks of 
oscillation in L for a fixed E increases rapidly as we go 
to lower values of L. When this distance is very small 
compared to the resolution width of L, the oscillation 
effect is averaged out. Only when the distance is large, 
it contributes to oscillation measurements. To get the 
reflection of this oscillation pattern in % 2 , we need de- 
creasing bin size of L with decreasing of its value. This 
has been studied in detail for three common choices of 
binning of the data in 15], and it has been found that 
neither log E — cos # ze nith nor in log E — log L is optimum. 
In this work we optimize the binning the data in the grids 
of log-E — L 0A plane. 

The number of bins used for this analysis is discussed 



later in the section IVll Here, it should be noted that one 
cannot make the bin size arbitrarily small. The number 
of event in a bin may be a fraction of 1 in theoretical 
data for chi-square analysis, but the number of event in 
experimental data is either zero or integer number greater 
or equal to 1. Obviously, no chi-squarc method will work 
if many of the bins have number of event equal to zero 
or just equal to 1. However, the number of events per 
bin > 1 is not also sufficient. We have checked that one 
needs number of event per bin at least > 4 to obtain 
% 2 /d.o.fw 1. This indicates the optimization of bin size 
with statistics. 



V. SELECTION OF EVENTS 

The up going and down going events are mixed at the 
near horizon due the uncertainty in scattering angle be- 
tween neutrino and muon. The up going neutrinos get 
oscillated and down going neutrinos remain almost un- 
oscillated due to the short distance from the source to 
the detector. When the iron plates of ICAL detector are 
placed horizontally, all these events cannot be detected. 
The high energy events will have normally small scat- 
tering angle, but very long tracks in the detector. So, 
they may be detected. If we plot the distribution of the 
difference in zenith angles between neutrino and the cor- 
responding muon for a fixed energy, it gives a Gaussian 
plot. The half width at the half maxima of this distri- 
bution as a function of muon energy is shown in fig [3] 
We put a selection criteria that the events for a given 
muon energy having the difference |90° — # ZC nith| within 
the above half width are rejected. Here we expect roughly 
that these events cannot be detected in case of real ex- 
periments. The precisions determined with and without 
this cut is discussed later. 



VI. THE x 2 

The number of events falls very rapidly with increase 
in energy and there is a very small statistics at the high 
energy. However, the contributions to the sensitivities of 
the oscillation parameters is significant from these high 
energy events. For the low statistics at the high energy, 
the x 2 is calculated according to the Poisson probability 
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FIG. 2: The — > oscillation probability in vacuum. We choose Am| 2 = -2.5 x 10 3 eV 2 , 6» 2 3 = 45° and 6\ 3 = 0° 
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FIG. 3: The variation of the half width at half maxima with for the distribution of #„ zcnith - # Al zcnith at horizon. The distribution 
is obtained for each i? M bin from 500 years un-oscillated atmospheric data of 1 Mton ICAL. 



distribution defined by the expression: 



X 



E 



2 K 1+ E^ 



k=l 



-2N° In 



(7) 



Here, N°j is the number of observed events gener- 
ated by Nuance for a given set of oscillation parameters 
with an exposure of 1 Mton. year of ICAL and Nfj is 
the number of predicted events (discussed later). These 
are obtained in a 2-dimensional grids in the plane of 
\ogE — L 0A . The term is the systematic uncertainty 
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of Nf } due to the fcth uncertainty (discussed later) and 
£k is the pull variable for the fcth systematic uncertainty. 
We use total number of log 12 bins n E = 35 (0.8 — 40 
GeV) and the number of L 0A bins as a function of the 
energy. We consider n L = 2 x 25, 2 x 27, 2 x 29, 2 x 31, 
and 2 X 33 for E = 0.8 - 1, 1-2, 2-3, 3-4, and > 4 
GeV, respectively. For the down-going events, the bin- 
ning is done by replacing 'L 4 ' by ' — L 0A \ The factor 
'2' is taken to consider both up and down going cases. 
For the up going neutrino, L is the distance traveled by 
the neutrino from the source at the atmosphere to the 
detector in the underground. In case of the down going 
neutrino, the L is the 'mirror L' which is the same L if 
the neutrino comes from exactly opposite direction. 

The table of Honda flux is given in 20 cos # zen ith bins 
and 101 E bins (0.1 - 10 4 GeV). It should be noted that 
we first re-binned the data into 300 cos Zenith bins and 
200 log E bins (0.8-40 GeV) to get the oscillation pattern 
accurately. This large number of cos # ZGn ith bins also help 
in proper re-binning of the data into L aA bins. 



A. Migration from neutrino to muon 



We have done this study for ideal muon detector. From 
GEANT simulation of ICAL detector it is seen that the 
energy resolution of muon varies 4-10% depending on the 
direction and energy. Since the iron plates are stacked 
horizontally, the resolution will be better for vertical 
events than the slanted events. The angular resolution 
varies from 4-12% for the considered range of energy and 
zenith angle. Here, the thickness of iron plates are con- 
sidered to be 6 cm. From Fig. [4] it is clear that these 
are negligible compared to the resolutions obtained from 
kinematics of scattering processes. 

The addition of the hadron energy to the muon energy 
of an event, which might improve the reconstructed neu- 
trino energy resolution, is not considered here for conser- 
vative estimation of the sensitivity. It would be realistic 
in case of GEANT-based studies since the number of hits 
produced by the hadron shower strongly depends on the 
thickness of iron layers. However, ICAL can also detect 
the neutral current events. Though it is expected that 
these events will not have any directional information; 
energy dependency of the oscillation, averaged over all 
directions can also contribute to the total x 2 m the sen- 
sitivity studies separately. 



To generate the theoretical data for the chi-square 
analysis, we first generate 500 years un-oscillated data 
for 1 Mton detector by Nuance. From this data we find 
the energy-angle correlated resolutions (see Figs. [4]) in 35 
E v bins (in log scale for the range of 0.8 — 40 GeV) and 
17 cos#^ enith bins (for the range —1 to +1). For a given 
\o%E v bin, we calculate the efficiency of having E^ > 0.8 
GeV (threshold of the detector). For each set of oscilla- 
tion parameters, we integrate the oscillated atmospheric 
neutrino flux folding the cross section, the exposure time, 
the target mass, the efficiency and the resolution function 
to obtain the predicted data. We use the CC cross section 
of Nuance-v3 [3] and the Honda flux in 3-dimensional 
scheme 
in 



B. Systematic uncertainties 

The atmospheric neutrino flux is not known precisely, 
there are huge uncertainties in its estimation. We may 
divide them into two categories: I) overall uncertain- 
ties (which are independent of energy and zenith angle), 
and II) tilt uncertainties (which are dependent of energy 
and/or zenith angle). We consider the following types of 
uncertainties. 

The energy dependent uncertainty, which arises due to 
the uncertainty in spectral indices, can be expressed as 



18| . This method has been discussed in detail $ Se (E) — 3>q(E) 



14j , but the number of bins and resolution functions 



EoJ 



ME) 



1 + S E log 



10 



E 
Eq 



have been changed here. 

One can do this directly by generating 500 Mton. year 
data (to ensure that the statistical error is negligible) 
for each set of oscillation parameters and then reduc- 
ing it to lMton.year equivalent data, which would be 
the more straight forward method. The marginalization 
study with this method is almost an un-doable job in a 
normal CPU. However, an exactly equivalent result is ob- 
tained here using the energy-angle correlated resolution 
function. 



(8) 

Similarly, the vertical/horizontal flux uncertainty as a 
function of zenith angle can be expressed as 

3«„(cob0,) « $o(cosfl 2 ) [l + tf.fl cos6 z \ - 0.5)] . (9) 

Next, we consider the overall flux normalization uncer- 
tainty 6 f N , and the overall neutrino cross section uncer- 
tainty 8 a . 

For E < 1 GeV we consider Se — 5% and Eq = 1 
GeV and for E > 10 GeV, S E = 5% and E = 10 GeV. 
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We take S fN = 10%, S a = 15%. We consider S z = 4% 
which leads to 2% vertical/horizontal flux uncertainty. 
We derived these uncertainties from [24J . 

For each set of oscillation parameters, we calculate the 
X 2 in two stages. First we used such tha t 4^ = 0, 
which can be obtained solving the equations [271 ] . Then 
we calculate the final \ 2 with these values. Finally, we 
minimize the \ 2 with respect to all oscillation parameters 
0. 

VII. MARGIN ALIZATION AND RESULTS 

A global scan of x 2 is carried out over the oscillation 
parameters Am 2 2 , 2 3, #13 and Sqp with neutrinos and 
anti-neutrinos separately. We have chosen the range of 
|Ato 2 2 | = 2.0-3.0xlCT 3 eV 2 , 23 = 38° -52°, 13 = 0°- 
12.5° and 8 C p = 0°-360°. The 2-dimensional 68%, 90%, 
99% confidence level allowed parameter spaces (APSs) 
are obtained by considering \ 2 = Xmin + 2-48, 4.83, 9.43, 
respectively. For every set of data we have checked that 
chi-square/d.o.f remains 1.1 at its minimum value. We 
obtain the APS in |Am| 2 | — 6*23 and |Amf 2 | — 6*13 plane. 
We set the input of |Am 2 2 | = 2.5x 10~ 3 eV 2 and 5 C p = 0. 

It is important to note here that the statistics changes 
significantly over L — E plane with the change of oscilla- 
tion parameters. Moreover, the fluxes and the resolutions 
are very different at different L — E zones. The upper and 
lower bounds of an oscillation parameter depends signif- 
icantly on the statistics as well as on the resolutions of 
the specific zones in L — E plane. The binning of the 
data, which captures the oscillation patten also plays the 
vital role. 

However, for some sets of input parameters the chi- 
square remains almost flat over a significant range of a 
parameter and then changes rapidly. It happens due to 
the fact that I) the change of oscillation probability is 
insignificant, and/or II) the above change is significant, 
but it is eaten by the systematic uncertainties in chi- 
square analysis. In this circumstances, the best-fit values 
may change significantly from the input values. This 
is a very common feature in analyses with generating 
events by Monte Carlo method. But, in methods without 
Monte Carlo, the number of events are determined with 
an accuracy of a fraction of 1 and then best-fit values is 
always close to the input values. 

In some cases, the deviations of the best-fit values 
are large. This is happened due to the following rea- 
sons. Here, we have just folded the total charge current 



cross section of all processes to find the number events 
for a particular neutrino energy to generate the theoret- 
ical data. We see significant fluctuations more than 1 
a in number of events between "theoretical data" and 
"experimental data" in some particular energy bins for 
a given set of oscillation parameters (see Fig. [5J. This 
happens mainly at the neutrino energy £ 3 GeV, where 
the resonances occur. Here, the neutrino cross sections 
depend on the type of nucleus. The generation of events 
is very complicated here and it depends on the mod- 
els. These all are not considered in the same way as in 
Nuance in generation of theoretical data, which causes 
energy dependent systematic uncertainty. However, this 
has no regular pattern. In our analysis we consider only 
the over all uncertainty in the cross section. These en- 
ergy dependent uncertainties have not been considered 
in our analysis. When 2 3 deviates from 7r/4, the differ- 
ence between peak and dip decreases and the fluctuations 
becomes relatively prominent. Again, when #13 becomes 
large, the periodic pattern of oscillation is lost due to 
matter effect. We have checked that the fluctuations are 
larger for 23 = 50° and 9 13 = 7.5° than 23 = 45° and 
#13 = 0°. In this region of oscillation parameters, signifi- 
cant deviations of best-fit values of oscillation parameters 
from their true values are obtained. 

The variation of A\ 2 [= X 2 ~ Xminl w ^h each of 623 , #13 
and Ara 2 2 , are shown in Fig. El El and [9] These are 
after marginalization over all the oscillation parameters 
except one, with which it varies. We present the cases 
for inputs 0i3 = 0°,and 7.5° with 23 = 40° ,45°, and 
50°, respectively. 

A. Sensitivity to #23 and its octant discrimination 

As the present experiments indicate that the value of 
013 is small compared to 023, the atmospheric neutrino 
oscillation is mainly governed by two flavor oscillation 
v v iPii) v t ipr)- This constrains sin 2 20 2 3 and |Am 2 2 |. 

From Fig. [TQl we see that the deviation from the max- 
imal mixing between 2 and 3 flavor eigen states can be 
observed. However, a degeneracy in 2 3 arises in case of 
013 = 0, whether it is larger or smaller than 45°. But, 
when the matter effect comes into the play, a resonance 
occurs in i/„ (p^) <-> v e (p e ) oscillation and it leads to 
a large effective value of 0i3 (see Eq. [4]). This helps to 
dominate the sin 4 2 3 term in Eq. [3] and breaks the 023 
degeneracy in its measurement. Since the atmospheric 
neutrinos cover a large region of E — L plane, it can ob- 
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FIG. 4: The sample energy-angle correlated resolution plots for neutrino (left column) and anti-neutrino (right column) for the bins 
of E v = 0.85 - 0.98 GeV with cos0 zc „ it h = -0.40 to -0.20 (upper row) and E v = 6.84 - 7.86 GeV with cos6» zoni th = to 0.20 
(lower row). The data is obtained from the simulation of 500 MTon.year exposure of ICAL considering no oscillation. 



serve the matter resonance and has an ability to discrim- 
inate the octant degeneracy. In Fig. [HI the variations of 
Ax 2 — (x 2 — Xmin) with #23 are shown for input values 
of 6> 23 = 40°, 45° and 50° with 6 13 = 0° and 7.5°, re- 
spectively. We see that with increase in 13 , the matter 
effect not only discriminates the octant, but increases the 
precision also. 



In Fig. [10] we see that for 613 = 7.5° the octant dis- 
crimination is possible for input of #23 = 40° and 50° 
with IH. But it is not possible for NH. Normally, the flux 
of Ufj, is higher than D^. In case of IH (NH), (2/^) is sup- 
pressed. The statistics remains high for IH compared to 
NH, which leads better octant discrimination possibility 
for IH. 



B. Sensitivity to #13 

The effect of #13 in oscillation probability does not ap- 
pear dominantly neither in atmospheric nor in solar neu- 
trino oscillation, but as a subleading in both oscillations. 
In case of atmospheric neutrino, its effect is seen at a) 
E ~ 1 GeV for propagation of neutrinos through vac- 
uum as well as through matter (no matter resonance), 
and b) E fts 2 — 10 GeV for propagation only through 
matter (matter resonance). The matter effect enhances 
the difference in oscillation probabilities between two #i 3 
values for neutrinos with NH and for anti-neutrinos with 
IH (see Eq. [4|. In Fig. [7] we show the cases a) and b) 
considering neutrinos and anti-neutrinos separately. We 
find that the effect of case a) is negligible. 

We have plotted the APS in #i 3 — |Am| 2 | plane in Fig 
rrrjfor e 13 = 0°, 5°, and 7.5° with 6> 23 = 40°, 45° and 50°, 
respectively. We find that the matter effect significantly 
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Neutrino (Experiment) — 
Neutrino (Theory) 

Anti-neutrino (Experiment) 

Anti-neutrino (Theory) 
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cos9 zenith 



FIG. 5: The typical distribution of events with E v keeping cos Zenith fixed at m —0.367 and with cos Zenith keeping E v fixed at 
« 2.24 GeV. We set Am§ 2 = -2.5 x KT 3 eV 2 , 6 23 = 45°, 13 = 0° and 5 C p = 0°. 



constrains #13 over the present limit. Though the matter 
effect acts either on neutrinos or on anti-neutrinos de- 
pending on the type of the hierarchy, but we have checked 
that it improves when we consider both neutrinos and 
anti-neutrinos. The sensitivity of #13 is not generally 
expected to be improved for the case of analysis with 
neutrinos and anti-neutrinos in together. However, this 
happens here due to the marginalization which restricts 
#23 more tightly for the case of v and v in together than 
either with v or D and indirectly constrains #13. It is also 
seen that the APS is strongly dependent on the input 
of 023 and a better constraint is obtained for #23 > 45°. 
However, it is notable here that the uncertainty is very 
high and the best-fit values deviate largely from its input 
values for nonzero 6*13 inputs due to the reasons discussed 



at the beginning of this section. 

C. Sensitivity to Am^ 

We show the constraint on |Am3 2 | in Fig. fTTJI and ITTI 
We see that the precision is little better when #13 = 0. 
The reason behind this is that a regular oscillation pat- 
tern with periodic rise and fall is observed when #13 = 0. 

It is seen that the APS is larger for NH than IH. The 
matter effect does not act on neutrino for IH and anti- 
neutrinos for NH with an addition to the fact that the flux 
is higher for neutrino than anti-neutrino. As discussed 
above, the APS is more restricted when there is no matter 
effect. Here, for input with IH the number of neutrino 
events is high and they do not have any matter effect. 
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Input IH, e 13 =0° Input 9 13 =7.5° 




respectively. Here, we consider both neutrinos and anti-neutrinos together. The type of input hierarchy is inverted. 



Input 9 13 =0°, 9 23 =45° 
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FIG. 7: The variation of A\ 2 = (x 2 ~ Xmin) w 'th $13 f° r input value of 9i3 = 0°, and 623 = 45° considering neutrino, anti-neutrino 
and both types of neutrinos, respectively. 



This leads to smaller APS for IH compared to NH fol- 
iar ge values of ^13. 

D. Effect of events at near horizon on precision 
measurements 

For a given set of input parameters, if we compare the 
APSs with zenith angle cut (discussed in sectionfVj with 
those without any cut, we find no significant differences. 
As a demonstrating example, we have shown the APSs 
in Fig. [12] without imposing any zenith angle cut for 
a given set of oscillation parameters. One can find the 
corresponding plots with zenith angle cut in Figs. 1 101 and 

EH 

From the study of this paper, we can conclude that the 
events at near horizon cannot contribute significantly in 
precision measurements. The fact is that the L resolution 
is very poor here. A little change in zenith angle at near 
horizon changes L values drastically. Again, the discrim- 
ination of up and down going events are not possible. 
So, the oscillation effect is almost smeared out by the 
resolutions. From the L/E dip considering the L and E 
values of neutrinos, one can expect a large contribution 



in precision from these events. But, in practical situa- 
tion, there is no appreciable improvement after addition 
of these events. 

The vertical (horizontal) stacking of iron plates will be 
able to detect the horizontal (vertical) events. So, from 
this study one can conclude that horizontal stacking is 
expected to give better precision than the vertical stack- 
ing. 

E. Precision of the parameters 

For a quantitative assessment of the result, we define 
the precision of a parameter t as: 

/ t ma,x_ t min\ 
P = 2[ r- . 10 

We find that the precisions are strongly dependent on 
the set of input parameters. We obtained the precision of 
|Am 2 2 | « 6.4%, 8.8% and 12% at 68%, 90% and 99% CL, 
respectively and the precision of sin 2 #23 ~ 31%, 38%, 
and 41% at 68%, 90% and 99% CL, respectively for the 
input of #23 = 45° and 9 13 = 0°. 
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FIG. 8: The same as Fig. [6] but with 6 13 . 



Input e 


13=0° 










B 23 
fl? 3 . 


=40° /— 
=4£°/ 

=ssn/ 





































Input 1 







e 23 =4oW— — 














1 





































i- CM 



FIG. 9: The same as Fig. [6] but with A 
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The oscillation dip moves towards the lower L/E val- 
ues as |Am| 2 | increases. The statistics also decreases at 
the lower L /E region. So, the precision is expected to be 
weaker as the input of |Amg 2 | increases. 

A comparison of the precisions of Arn| 2 an d sin 2 # 23 
among different future baseline experiments is made in 
[3]. The variation of the precisions with the change of 
input parameters are also presented there. We compare 
our results with 5 years run of T2K, which is the best in 
determining precision of atmospheric oscillation param- 
eters in the list in [3]. The precision of Am§ 2 is almost 
same with T2K (w 12%) and precision of sin 2 823 from 
ICAL is 41% while from T2K is 46%. Here we present 
the results for 10 years run of 100 kTon ICAL detector. 
From this work it is also seen that atmospheric neutrinos 
at ICAL detector are in very good position to discrim- 
inate octant of 623- The main advantage here is that 
atmospheric neutrinos are natural sources and the cost 
goes only to build and run the detector. 



VIII. CONCLUSION 

We have studied the precisions of the oscillation pa- 
rameters from atmospheric neutrino oscillation experi- 



ment at the large magnetized ICAL detector generating 
events by Nuance and considering only the muons pro- 
duced by the charge current interactions. The distance 
between two consecutive peaks of oscillation in E for fixed 
L increases as one goes from higher L values to its lower 
values. This indicates the need of finer binning at lower 
L values in \ 2 analysis. We optimize the binning of the 
data in the grids of \ogE — L QA plane. We find that 
the impact of the events at near horizon on the precision 
measurements is very negligible due to poor L resolution. 

From the marginalized \ 2 study separately for neu- 
trinos and anti-neutrinos, we find that the measurement 
of #13 is possible at a considerable precision with atmo- 
spheric neutrinos. The precision of 613 depends crucially 
on its input value. For 613 = 0, we find its upper bound 
ps 4°, 6° and 9° at 68%, 90% and 99% CL, respectively. 
The both lower and upper bounds of 613 are also pos- 
sible for some combinations of (#23,^13) and it happens 
mainly for 623 ^45°. 

The precision of |Am§ 2 | and 823 can also be very high 
and the determination of octant of 823 is possible for some 
combinations of (#23, #13 )■ 

It should also be noted here that in x 2 analysis the 
theoretical data and the experimental data are not gen- 
erated in the same way. The different models of neutrino 
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interactions generate energy dependent systematic un- 
certainties at some energies. These are not included in 
this analysis. This causes sometimes large deviation of 
the best-fit values of the oscillation parameters from the 
input values. 
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FIG. 10: The 68%, 90%, 99% CL allowed regions in #23 - |Am| 2 | plane for the input of 6» 23 = 40° (first row), 45° (second row), 
50° (third row) with 6*13 = 0° (first column), 5° (second column), 7.5° (third column) with IH and 7.5° (fourth column) with NH. 






FIG. 11: The same as Fig. \T6\ but in 9i3 - |Ara| 2 | plane. 




FIG. 12: The allowed regions without any zenith angle cut for the events at the horizon. 



